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Abstract 

This study of spatio-temporal pattern formation in an annulus is motivated by 
two physical problems on vastly different scales. The first is atmospheric convec- 
tion in the equatorial plane between the warm surface of the Earth and the cold 
tropopause, modeled by the two dimensional Boussinesq equations. The second is 
annular electroconvection in a thin smectic film, where experiments reveal the birth 
of convection-like vortices in the plane as the electric field intensity is increased. This 
is modeled by two dimensional Navier-Stokes equations coupled with a simplified 
version of Maxwell's equations. The two models share fundamental mathematical 
properties and satisfy the prerequisites for application of 0(2)-equivariant bifurca- 
tion theory. We show this can give predictions of interesting dynamics, including 
stationary and spatio-temporal patterns. 

Keywords: Annular Rayleigh-Benard convection; Annular electroconvection; 
Bifurcation; Steady-state mode interactios; 0(2) symmetry. 



1 Introduction 

Pattern formation is observed in a wide variety of contexts, from slime 
moulds to fluid mechanics [1,2]. There is a growing realization that certain 
aspects of pattern formation phenomena are model independent, that is to 
say, independent of physical details of the particular model. Instead, the types 
of patterns which may be observed are largely determined by purely mathe- 
matical properties such as geometry and symmetry. As an illustration of this 
principle, we consider two different examples of pattern formation, which differ 
physically in the source of the instability (one thermal and the other electrical) 
and occur on vastly different spatial scales (ratio ~ 10^). What the two exam- 
ples have in common is their geometry and symmetry groups. Both take place 
in an annular domain and have the symmetries of the group 0(2). This paper 
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presents an analysis of those aspects of pattern formation in the two examples 
which are model independent. In the main, these results are qualitative, not 
quantitative. 

Of course, a more detailed description of the behavior of a model is possible 
if one takes into account all of its physical properties. The extension of the 
qualitative analysis presented here to a detailed quantitative analysis is in 
progress and will be presented elsewhere. This type of quantitative analysis 
was first carried out in the context of the Taylor-Couette problem [3,4], which 
shares much with the present examples in terms of geometry and symmetry. 
Other examples where 0(2) symmetry plays an essential role include buckling 
of a cylinder, fiow through a cylindrical pipe and systems which are periodic on 
a line. The separation of quantitative from qualitative analysis is intended to 
illuminate the distinction between model independent phenomena and those 
which are intrinsically physical. The mathematical tools used here, known as 
Equivariant Bifurcation Theory, have been developed over the past twenty 
years by many mathematicians and physicists; good references are [1,2,4]. 

In brief outline, this paper presents two different examples of annular con- 
vection and extracts their common features on which a model independent 
analysis may proceed. A summary of equivariant bifurcation theory is pre- 
sented, with emphasis on the key features of equivariant mode interactions in 
the case of 0(2) symmetry. The predictions of this mathematical theory for 
pattern formation in annular convection are presented in tables and bifurca- 
tion diagrams. These possibilities include stationary and spatio-temporal pat- 
terns, hysteresis and mode-jumping, smooth symmetry-breaking transitions, 
rotating waves and standing waves. 



2 Annular Convection: Two Problems 

An annulus is the region between two concentric circles in a plane. We 
consider an annulus filled with an incompressible fiuid, subjected to an applied 
force vector directed radially with magnitude independent of the angle, see 
Figure 1(a). If this radial force is sufficiently intense, the fiuid may begin to 
move in interesting patterns as in Figure 1(b) and (c). A fundamental problem 
is how to predict these patterns. 

An objection which may be raised on physical grounds is that this annulus 
is two dimensional, and physically no fiuid is two dimensional. However, in 
both of the examples presented here, there are strong arguments for a two 
dimensional model of a three dimensional fiuid. Special features of both ex- 
amples, in the one case geometrical and in the other case physical, cause the 
fiuid to behave as if it were two dimensional. This greatly reduces the com- 



2 



putational difficulties in the models, and facilitates the analysis of pattern 
formation. Also, they provide a setting in which predictions from a two di- 
mensional mathematical fluid model can be compared with physical reality. 

2.1 Thermal Convection in the Equatorial Plane 

Imagine that the Earth's atmosphere is spherical and ignore imperfection ef- 
fects, such as continents and the seasonal tilt, which break the symmetry of 
a sphere rotating about a North-South axis. Then the plane slicing the Earth 
at the equator is a plane of symmetry and there are no forces acting on the 
atmosphere normal to this plane, only tangential forces. This idealized sym- 
metry implies the possibility of flow patterns for which the equatorial plane is 
invariant; that is, there is no flow across the equatorial plane. Eor the actual 
Earth, it is known that this is very nearly the case. 

All clouds and other weather phenomena of the atmosphere are to be found 
in the troposphere, which is bounded below by the warm surface of the Earth 
and bounded above by the cold tropopause. The equatorial plane intersects 
the troposphere in an annular domain. We restrict our attention to the atmo- 
spheric flow in this annulus around the equator. Recent observational work in 
geophysical fluid dynamics has revealed unexplained large-scale waves in the 
equatorial troposphere [5,6]. 

The geophysical fluid dynamics (GED) in this atmospheric annulus may be 
modelled accurately by the Boussinesq system of equations [7]. On a thermal 
diffusion time scale, the system is written in the nondimensionalized form 
(after a translation T = T — T^^ p = p — p^^ dropping the ~ and neglecting 
the Coriolis and centrifugal forces due to the rotation of the Earth) 



where the variables are 
t time, 

(r, 6) polar coordinates, 

V nabla operator in polar coordinates: V = (d/dr^ r~^ d/dO)^ 

V fluid velocity: v = {vy{r,6 ,t)^ve{r,6 ,t))^ 



dv 
'dt 



+ (iT . V) tT= Pr [-Vp + V^v + Ra (Tg)] , 
V • u = 0. 




dT 
'dt 
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T temperature variation: T = T{r,6,t), 
Tg temperature of equilibrium: Tg = Tf.{r), 
g gravity unit vector in the radial direction, 
Pr Prandtl number, 
Ra thermal Rayleigh number. 

The boundary conditions on the two circles are 

= 0, r = 0. 

2.2 Electroconvection in a Smectic Liquid Crystal Film 

There exist thin liquid crystal films (of type smectic- yl) which behave like an 
isotropic incompressible liquid in a plane but like a plastic crystal in the normal 
direction [8] . The film is formed of a small integer number of planar molecular 
layers, of uniform thickness and behavior throughout. This structure fixes the 
plane of the smectic layers and suppresses any fiow normal to the film. The 
film also is a weakly conducting dielectric, exhibiting an ohmic behavior for 
low electric field intensities in the plane. Experiments reported by S. Morris 
[9] on a thin smectic- yl liquid crystal film suspended in a small annulus (order 
of 10~^m), with an electric field applied in the radial direction, reveal the birth 
of convection-like adjacent pairs of counter-rotating vortices in the plane as 
the electric field intensity is increased, see Figure 1. 

In [10] a surface charge model is developed for these electrohydrodynamic 
(EHD) experiments, the instability being interpreted in terms of the inter- 
action between the exterior three dimensional electric field (just outside the 
film) and the surface charges developed on the smectic- yl film viewed as a 
suspended surface with two (top and bottom) faces. In this model, it is the 
discontinuity in the normal component of the electric field across the surfaces 
of the smectic film which creates surface charges. These surface charges in 
turn experience a force due to the tangential component of the electric field, 
which drives the convection instability. 

On the other hand, a volume charge model may be constructed, taking into 
account the spontaneous generation of ions in the bulk of the liquid crystal 
through electrochemical dissociation/recombination, perhaps involving dopant, 
(see [11]) 

AB ^ A- + B+. 
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This dissociation by itself yields equal numbers of positive and negative ions 
and a zero net charge; however, near the positive and negative electrodes re- 
spectively positive and negative net charges build up, due to injection and/or 
annihilation mechanisms. Locally (near the electrodes) diffusion is important 
in this model. A convective instability arises from the interaction between 
these net volume charges and the two-dimensional electric field tangent to 
the film, analogous to that in the Boussinesq model. Experiments reported in 
[8] on rectangular smectic-yl films, with and without electrode contact, sug- 
gest that some contact mechanism at the electrodes, perhaps electrochemical, 
is required for pattern formation. The following equations represent a vol- 
ume charge model which neglects the surface charge mechanism, and strongly 
resembles the above Boussinesq approximation in GFD. After scaling on a 
conduction charge decay time, it may be written in the nondimensionalized 

form (after a translation q = q — q^^ E = E — E^^ p = p — and dropping 
the ~ ) as 

-^ + (iJ.V)v = PT [-Vp + V^v + Rci [qE + qJ + qE,)] , (2) 
V • u = 0, 

Oq 

— + (v ■ V) q = q + S Pt V'^q - ql{r) v^, 

q = VE, 
E = -V<t>, 

where the variables are 

q net charge density variation: q = q{r,6,t), 

q^ net charge density at equilibrium: q^ = qe{r)^ 

E in-plane electric field variation: E = {Ey{r, 6, t), Ee{r, 6, t), 

Ef. in-plane electric field of equilibrium: E^. = (£'e(?"),0), 

cf) electric potential, 

Pr electric Prandtl number, 

Ra electric Rayleigh number, 

8 dimensionless ratio of charge diffusion coefficient to kinematic viscosity. 

The boundary conditions on the two circles are 
u = 0, q = (f> = ^. 
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If two or more species of ions with different mobilities are important, then this 
model may be expanded with a different charge continuity equation for each 
species of ion. In the forcing term of the Navier-Stokes equation the overall 
total net charge will appear. 

2.3 Common Features of the two Models 

In spite of their obvious differences, the two annular convection models from 
GFD (geophysical fluid dynamics) and EHD (electrohydrodynamics) have 
many common features. Both are especially interesting because they provide 
physically relevant examples of two dimensional Navier-Stokes flows. Both 
have 0(2)-symmetry, i.e. the symmetry of the circle under rotations and re- 
flections across any diameter. 

There is a destabilizing radial force, due to gravity and the temperature- 
induced density variation for the GFD model and due to the in-plane electric 
field intensity and the electric charge density variation for the EHD model. 
Each of these destabilizing forces is opposed by a a stabilizing force due to 
viscosity. 

In both cases there are trivial (pre-convection) solutions which are known 
exactly (analytically). The neutral stability boundary is fairly easy to calculate 
(numerically) after an exact linearization about the known trivial solution. 
The initial loss of stability (bifurcation) is to a convection mode, which can 
be described as a chain of pairs of clockwise/counterclockwise vortices strung 
periodically around the annulus, see Figure I. The number of such pairs is 
the mode number m, and this depends on the radius ratio rj G (0, 1). The two 
important parameters are the Rayleigh number Ra (depending on temperature 
or voltage respectively, as appropriate to the model) and the radius ratio rj. 
Codimension-two points for two different mode numbers exist in the parameter 
space. We show that the prerequisites for the application of 0(2)-Equivariant 
Bifurcation Theory are satisfied for both models. This implies the possibility 
of rich mode interactions. 



3 Equivariant Bifurcation Theory 

We investigate these two problems in parallel, since the mathematical 
tools [f-4] for both are essentially the same. Consider a parametrized family 
of differential equations, representing for example (I) or (2), 

y' = f{y,^i), f:XxW^y, (3) 
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where / is a smooth function {k > 2) such that /(0,0) =0, X,y are 
appropriate function spaces and TV is the parameter space [TZ = set of real 
numbers). Local bifurcation theory deals with the question: What qualita- 
tive changes can happen to the set of solutions near j/ = 0, as the parameter 
/i varies near zero? The answer depends in large part on the eigenvalues of 
the linearization of / at the origin, i.e. Dy/(0,/i). Typically there are three 
possibilities: either there is no qualitative change (the eigenvalues remain hy- 
perbolic as /i varies), or there is a steady-state bifurcation (a real eigenvalue 
satisfies A(/i) = for some //), or there is a Hopf bifurcation of a limit cycle 
(a complex eigenvalue satisfies Re [A(/i)] = for some //). In this paper, we 
consider only the case of steady-state bifurcation, but we allow the possibility 
of multiple real eigenvalues vanishing for the same value of i.e. of mode- 
interactions. Without loss of generality, we translate fi so that this multiple 
zero eigenvalue occurs at // = 0, i.e. A(0) = 0. 

3.1 Reduction to a Center Manifold 

Suppose that this zero eigenvalue A(0) = has finite multiplicity and that its 
complement in the spectrum of Dy/(0, 0) is bounded away from the imagi- 
nary axis, in the negative half-plane. Then, under reasonable conditions on /, 
which are satisfied for many PDE of the general form (3), there exists a center 
manifold in A! which is tangent to, and of the same dimension as, the (general- 
ized) eigenspace of this zero eigenvalue [12]. Even in cases where the existence 
of this center manifold has not yet been proven rigorously, it is possible to 
calculate it in a formal manner. In this paper, we set aside the pure existence 
question and assume there is a normally hyperbolic center manifold for the 
annular convection models above. Then, locally, all solutions off the center 
manifold approach it exponentially with increasing time. Thus the study of 
the long-time behavior of solutions of (3) near the origin is reduced to that 
of solutions on this finite-dimensional manifold. This reduction to finite di- 
mensions has been calculated, for example, in the case of the Taylor-Couette 
problem in [3,4]. 

Furthermore, this reduction process preserves the symmetry properties of the 
system, as described in the next subsection. 

3.2 Symmetry: Definitions and Tools 

Consider a system of ordinary differential equations 

j/' = /(j/,/i), ^/G7^^ /iG7^^ (4) 
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and let F be a compact group acting linearly on 7^" i.e. 

^:7^"^7^", rcGL(n), 

where GL(n) is the group of all invertible nxn matrices over TZ. Then, system 

(4) is said to have the symmetry F (or is T-equivariant) if 

/(7?/,/^) = 7/(?/,/^), V7Gr, (^/,/i)G7^'^x7^^ (5) 

In this paper we are concerned with the orthogonal group T = 0(2) on TZ^ 
which is generated by Rg, the planar rotations by ^, and the flip /€, the re- 
flection across the Then, the definition (5) says that / has symmetry 
0(2) if / commutes with every 7 G 0(2). This algebraic identity agrees with 
our geometrical notion of symmetry. Indeed, taking 7 = Rg, the definition of 
F-equivariance says that the value of / at the point obtained by rotating y 
through angle 6 is the same as f{y) rotated through 6. Similarly, taking ^ = 
it says that / at the refiection of y is the refiection of f{y). 

For the remainder of this section, let us temporarily suppress the parameter 
/i in (4) and (5), since // plays no direct role in the symmetry; that is rewrite 

(5) as 

/(7x) = 7/(x), V7Gr, ^/G7^^ 

The group orbit of a point x G 7^", with respect to F, is defined as 

Tx =: {^x I 7 G F}. (6) 

Suppose that x is a steady-state, i.e. f{x) = 0. Then = 7O = ^f{x) = f[^fx)^ 
for all 7 G F, so the entire group orbit is composed of steady-states. We have 
the following dichotomy (and both cases are interesting): either ^fx = x (i.e. 
X has the symmetry F) oi ^fx ^ x (i.e. 7a; is a new solution). In the latter 
case, each 7a; is called a conjugate of x. Conjugates have the same dynamic 
behavior. Similar considerations may be employed for limit cycles and other 
types of solutions of (4). 

The isotropy subgroup of x G TZ^ is defined as 

Sx =: {7 G F I -fx = x} (7) 

and represents the total symmetry of x. The fixed-point subspace of E, defined 
by 

Fix(E) =: {x G 7^" I 7x = X, V 7 G S} (8) 
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represents the set of all points fixed by the subgroup E. It is a subspace of 
7^", which is invariant for any F-equivariant /, i.e. 

/(Fix(E)) C Fix(E). (9) 

Even though / is nonlinear, we can restrict the problem to lower dimensional 
linear subspaces on which / is invariant. This is much like the decomposition 
of a linear problem into eigenspaces. Thus (9) yields a very powerful tool in 
the analysis of symmetric systems. 

Next we show how symmetry can lead to multiple eigenvalues. Assume that 
is a steady-state, i.e. /(O) = 0. Let L =: [df/dy){0) (the linearization of / at 
0), and suppose that L has the eigenvalue A = 0. Then V =: keiL ^ {0}. By 
differentiating the equivariance relation f{^fx) = ^f{x), we obtain, at x = 0, 
that = 7iv, so L and 7 are commuting linear operators. Thus F leaves V 
invariant. For most F, we must have dimV > I for a faithful representation 
on V . For instance, in order to have the full symmetry of F = 0(2) (the 
symmetry of the annulus) or F = Dn (the n-gon symmetry), we must have 
at least dimV = 2. Then A = is a multiple eigenvalue "forced" by the 
symmetry. 

3.3 Normal Forms 

For F compact, F-equivariant maps / can be written in a very useful simplified 
form. First consider the class of invariant functions, i.e. smooth functions 
g : TZ^ — > TZ satisfying 

g{^x) = g{x) V 7 G F, X en''. (10) 

Then there exists a finite set of invariant polynomials {ui, Us} with U]i{^fx) = 
U]i{x) for all 7 G F and k G {f, •••,5}, such that any invariant function g can 
be written as g{x) = h[ui[x), ...,Us{x)) for some smooth function h. This re- 
sult holds for g and h polynomials or germs. The set {mi,...,Ms} is called a 
Hilhert basis. Next follows the Equivariant Normal Form Theorem: for a given 
equivariant / : TZ^ — > (polynomial or germ) there exists a finite set of 
F-equivariant polynomial maps Vi, ...^v,- : — > such that 

r 

/(^) = h]^{ui{x), ...,Us{x))v]^{x) (II) 
k=i 

with /ik and Uj as above. Equation (II) gives the equivariant normal form 
of /. For the case of F = 0(2), the sets of basic invariant polynomials and 
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equivariants were obtained by Buzano and Russo, see [13,14]. 



3.4 Mode Interactions 

Roughly speaking, the codimension of a bifurcation is the number of equa- 
tions (in the parameters) which must be satisfied in order for that bifurcation 
to occur, under the assumption that certain nondegeneracy conditions hold. 
Typical bifurcations have codimension one. A simple setting in which a codi- 
mension two bifurcation can occur is as follows. In a two-parameter system, 
as we vary the two parameters fi = (/ii,/i2) G in y' = f[y^fi)^ we en- 
counter codimension-one curves of neutral stability, each corresponding to a 
zero eigenvalue (or complex eigenvalue with zero real part) with the symmetry 
of the corresponding kernel (see Figure 2). The curves corresponding to bifur- 
cations of codimension-one are well-understood. Generically these curves will 
intersect transversally at isolated points (codimension-two bifurcation points) 
where mode interactions occur. Near a codimension-two point as in Figure 
2 (at the intersection of two codimension-one curves) one naturally expects 
to find both of the codimension-one behaviors, and generally does so. How- 
ever, the interesting fact is that often there is much more to be found, since 
nonlinear "mode interactions" can produce interesting new phenomena. 



4 Analysis of Mode Interactions 

The analysis of the mode interactions takes place in the following steps. 
For more details, see [1]. 

• Adjust the parameters to a neighborhood of the codimension-two point of 
0(2) mode interactions, i.e where two eigenvalues are simultaneously zero. 
This involves only linear stability analysis. 

• Perform the center manifold reduction from the original PDF to a four 
dimensional ODF system 



(/(0,/i) = 0, D(/(0,0) = 0, (/(7j/,/i) = 7(/(j/,/i), V7G0(2) (13) 



(we assume 0(2) acts faithfully on both modes, which implies a four dimen- 
sional system). 






where 
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• Identify the state space TZ'^ with using the complex coordinates {zi,Z2). 
The action of 0(2) on is given by 9{zi,Z2) = (e^"^^ Zi,e^'^ Z2) (rotation) and 
n{zi,Z2) = (^1,^2) (reflection) where / and m are the integer mode numbers. 
Based on available experimental evidence, we may assume that / = m + 1. 

• The Hilbert basis for the ring of 0(2)-invariants as obtained in [13] is 
{u, given by 

u=:|zi|2, v=:\z2\^, w =: z^+^z^ + z^+^z:^ . (14) 



• The generators of the module of 0(2)-equivariants (over the above ring of 
invariants) from [13] is 



(15) 












■^m m 












Z2 


■} 
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m + 1— m-1 

^1 ^2 



• It follows that the normal form for the 0(2)-equivariant system is 



d 















^1 ^2 







dt 




= Pi 




+ P2 




+ qi 




+ 92 




Z2 





Z2 





m + 1— m-1 

Zl Z2 



(16) 



together with the two equations for Zi,Z2 which are complex conjugates of 
these. Here pi,P2, <?i, <?2 are real functions of (u, fi)^ and satisfy 

Pi(0) = P2(0) = 0. 

• Compute the isotropy subgroups of 0(2). They are conveniently organized 
in the form of a lattice of isotropy subgroups^ where the arrows represent in- 
clusions: 



0(2) 



D 



m+1 



Z2 

T 

{1} 



Table 1: Lattice of Isotropy Subgroups. 
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• In Table 1, Z2 is the reflection symmetry and {1} has no symmetry. Each 
isotropy subgroup determines a corresponding flxed-point subspace. 

• Restrict the system to each flxed-point subspace and solve it. We recall 
that each flxed-point subspace is invariant for the flow (i.e. for time-evolution 
of solutions) and the system simplifles when restricted to each flxed-point 
subspace, as follows: 



Isotropy 
Subgroup 
E 


Fixed-Point 
Subspace 
Fix(E) 


Equations Restricted to 
Eixed-Pomt Subspace 
5 Fix(S) = 


0(2) 




= 




{{x,0)\x G 7^} 


X pi{x^, 0, 0, /i) = 
zi = X €z TZ 

Z2 = 


Dm-l-1 


{(o,^/)|^/G7^} 


2/^2(0, ^/^0,/i) = 
zi = 
Z2 = y ^TZ 


Z2 


z = {x,y) 
x,y TZ 


hi=pi + qix'^-^y"' = 
h2=P2 + q2x"'+^y"'-^ = 


{1} 




Pi = P2 =0 

9i = 92 = 



Table 2: Restrictions to Fixed-Point Subspaces. 



• Compute numerical values of the leading coefflcients of (low degree) non- 
linear terms in , p2 , <?i , <?2 and thus determine local dynamics in each of the 
flxed-point subspaces in this table. (This is work in progress.) 



5 Bifurcation Scenarios 

A general analysis of the 0(2)-equivariant normal form (16) has been per- 
formed by Dangelmayr [14] and reviewed in [1]. A rich assortment of behaviors 
is possible. The determination of which behaviors actually occur in a given 
physical system depends on the speciflc values of certain coefflcients (to be 
computed). Here we describe some interesting possibilities which follow from 
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the general analysis, assuming m > 2. (The case m = 1 is quite different and 
will not be discussed here.) 



First, consider the equations for steady-states in the subspace Fix(Z2), i.e. 
fl'|Fix(Z2) = 0: 




0, 
0, 



(17) 



P2y + q2X^+^y^-^ = 



where pi, (j'i,p2, <?2 are functions oi u = ,v = , w = x'"^'^^y'"^ and the param- 
eters. It can be shown that, generically, there are no steady-states other than 
those in Fix(Z2). Thus symmetry has reduced the dimension of the problem 
from four to two. It is clear that (17) has the trivial solution (x^y) = (0,0). 
We seek nontrivial solutions. 

Suppose now that the mode number m is very large, as it would be in the GFD 
example (1) where the radius ratio is very close to 1. Then locally, the two 
terms including qi , q2 are very small compared to the other two terms involving 
P\,P2 (which are independent of m). Similarly, among the invariants, w = 
Re z'2 ) is very small for large m compared with u and v. Thus, to obtain 

a first approximation when m is large, we may set w^qi^q2 all equal to zero. 
The resulting equations are the well-known bifurcation equations for the case 
of Z2®Z2 symmetry, which arise in cases of pitchfork/pitchfork, pitchfork/Hopf 
and nonresonant Hopf/Hopf mode interactions (with no additional symmetry 
or degeneracy). The bifurcation equations for this case to leading order are 



where a, /3 are the unfolding parameters and pij are the obvious derivatives, 
evaluated at the origin. We assume the nondegeneracy conditions P11P22 7^ 
and P11P22 — P12P21 7^ 0. The complete set of bifurcation diagrams for this 
situation has been published in [15], see also [16,1]. Two interesting cases 
are shown in Figure 3, where the bifurcation parameter fi is related to the 
unfolding parameters by 

a = /i, l3 = fi -\- 6^ 

and (5 > is a small "splitting parameter". In all cases, there are two primary 
pitchfork branches from the trivial solution, given by 





y {(3 + P2ix^ + P22y^) = 0, 




(19) 
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x = 0, y = , 

P22 

and secondary "mixed-mode" branches with both of x, j/ ^ 0, which are solu- 
tions of the linear system in obtained from (18) on dividing out x ^ 
and y ^ 0. The existence of all these branches persists when the higher order 
terms are restored to the equations (18), by the Implicit Function Theorem. 
In the case of Z2 ® Z2 symmetry (with w^qi^q2 suppressed) the symmetries 
of these branches also are preserved: the primary branches remain symmet- 
ric pitchforks, and the secondary branches occur as four conjugate (iba;,ibj/) 
under the Z2 ® Z2 symmetry. It is convenient in this case to identify these 
conjugates and draw them as single branches in a planar bifurcation diagram 
with amplitude plotted against bifurcation parameter, see Figure 3. 

Now let us return to the equations for steady-states of the original 0(2)- 
symmetric system given by (17). The m-dependent terms break the Z2 ® Z2 
symmetry but the existence of all the branches in the Z2 ® Z2 symmetric case is 
preserved. For sufficiently large m, the symmetry-breaking effects will be very 
small and the bifurcation diagram will be indistinguishable from the Z2 ® Z2 
case. For all values of m > 2, we note that (17) retains a Z2 equivariance which 
depends on the parity of m: 

(x, J/) — )• (— x, J/), m odd 
(x, J/) — )• (x, — J/), m even . 

For definiteness, let us assume m is odd (the even case is similar). Then 
in the bifurcation diagrams of Figure 3 (full 3D versions), the Z2-symmetry: 
(x^y) — )• (—x^y) implies that the Dm-mode primary pitchfork is preserved 
and the secondary bifurcations occur at the same points on the two prongs 
of this pitchfork. However, the lack of a Z2-symmetry: (x^y) — )• (x, — j/) 
means that the four secondary branches occur in two distinct sets of conjugate 
pairs. Thus in drawing the conjugacy classes of branches in this case we show 
two distinct Z2-symmetric (mixed-mode) branches, as in Figure 4. A further 
implication of the breaking of the Z2-symmetry: (x, y) — )• (x, —y) is that the 
secondary bifurcation from the Dm-symmetric x-mode becomes an imperfect 
pitchfork, which implies the possibility of a transcritical bifurcation with a 
nearby saddlenode (or limit point) bifurcation. This is indicated in Figure 4. 
There is an additional "hidden symmetry" which preserves only the primary 
pitchfork for the Dm+i-symmetric j/-mode with x = 0, see [1] for details. 

The bifurcation diagrams in Figure 3 show secondary (mixed-mode) branches 
which link the two primary (pure-mode) branches, in two different ways. In 
the ), the secondary branch is itself unstable but it effects an exchange 

of stabilities between the primary branches (stable solutions are indicated by 
heavy lines). If the parameter fi is varied up and down through this range. 
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hysteretic jumps between the two primary branches will occur. In case (b) 
of Figure 3, the secondary branch is stable and it effects a smooth stable 
transition between the primary branches. Along the stable secondary branch 
joining the Dm to the Dm+i primary branch, the corresponding picture in the 
annulus has a small counter-rotating vortex pair created near the outer circle 
between two existing vortices. This new pair grows and eventually matches 
the other pairs as the Dm+i branch is approached. On this secondary branch 
the symmetry is broken, from (or Dm+i) to Z2. 

Additional phenomena can be predicted if one of the primary bifurcations is 
subcritical with the other remaining supercritical, as in Figure 4. Then, for 
the Z2 ® Z2 symmetric normal form, as shown in [16,15], there is generically 
a Hopf bifurcation of a periodic solution from the mixed-mode branch. In 
fact, there are four identical Hopf bifurcations but all are conjugate under the 
Z2 ® Z2 symmetry. Returning now to equations (17) for the steady-states of 
our original 0(2)-symmetric system, these Hopf bifurcations are preserved but 
now exist as two distinct pairs of Z2-conjugates. Thus, under these conditions 
we predict the Hopf bifurcation of two distinct branches of periodic solutions. 
Because these are purely temporal oscillations (not spatio-temporal) we may 
call them standing wave solutions. They are marked SW in Figure 4. 

Still another type of solution is possible, as a consequence of the original 0(2)- 
symmetry. The primary (pure-mode) and secondary (mixed-mode) solution 
branches found for the reduced system (17) actually exist as full "circles" 
of solutions, conjugate under the S0(2)-symmetry, for the original equations. 
Dangelmayr [14] pointed out that in addition to the two eigenvalues associated 
with the 2D-systems above, there is one eigenvalue identically zero due to this 
S0(2)-symmetry, and a fourth eigenvalue which may change sign. When this 
happens, the result is bifurcation of a rotating wave, involving a slow motion 
along the S0(2) group orbit. Here by rotating wave is meant a spatio-temporal 
pattern with the property that there is a rotating coordinate frame in which 
the pattern becomes stationary. This rotating wave is marked RW in Figure 
4. 



6 Conclusions 

A rich variety of stationary and spatio-temporal patterns is possible in 
systems with the symmetries of annular convection problems. In the geophys- 
ical fluid dynamics problem (1), the mode number m is very large, which 
implies that the bifurcation diagrams are practically indistinguishable from 
the classical Z2 ® Z2 symmetric case [16,15]. Fither hysteresis or smooth tran- 
sitions between modes, and stationary waves are possible behaviors. In the 
electroconvection problem (2), mode numbers m are in the single-digit range. 
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This allows even richer behaviors, for example the possibility of a rotating 
wave. Further calculations (in progress) will determine which of these can be 
predicted in given experiments. 
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FIGURE CAPTIONS 

Figure 1: Symmetry Breaking from Basic State, (a) Basic state, radial 
forcing, 0(2) symmetry; (b) Mode m = 3, three pairs of convection rolls; (c) 
Mode m = 4, four pairs of convection rolls. 

Figure 2: Neutral Stability Curves. Codimension-one curves intersect at 
codimension-two points in the (Ra, rj) plane, where Ra is Rayleigh number 
and T] is radius ratio. 

Figure 3: Sample Bifurcation Diagrams for Large Mode Number. 

Both primary bifurcations are assumed supercritical. The branches are labeled 
by their isotropy groups. Stable branches are indicated by heavy lines. 

Figure 4: Sample Bifurcation Diagrams for Small Mode Number. 

One primary branch is assumed subcritical and the other one supercritical. 
The labels SW and RW denote standing waves and rotating waves respectively 
and the other branches are labeled by their isotropy groups. 
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